Preheating - cosmic magnetic dynamo ? 
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We study the amplification of large-scale magnetic fields during preheating and inflation in several 
different models. Preheating can resonantly amplify seed fields on cosmological scales. In the 
presence of conductivity, however, the effect of resonance is typically weakened and the amplitude of 
produced magnetic fields depends sensitively on the evolution of conductivity during the preheating 
and thermalisation phases. In addition we discuss geometric magnetisation, where amplification 
of magnetic fields occurs through coupling to curvature invariants. This can be efficient during 
inflation due to a negative coupling instability. Finally we discuss the breaking of the conformal 
flatness of the background metric whereby magnetic fields can be stimulated through the growth of 
scalar metric perturbations during metric preheating. 
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I. INTRODUCTION 

With the current dominance of the inflationary 
paradigm and the gravitational instability picture of 
structure formation seeded by quantum fluctuations it 
is easy to forget earlier, competing, models. In partic- 
ular, models of structure formation based on turbulence 
had the advantage that they were able to make strong 
connections between galaxy formation, galactic angular 
momentum and galactic magnetic fields 

Inflation, by contrast, predicts essentially zero vortic- 
ity and in its purest forms^ rather small magnetic fields. 
The end of inflation may be very violent, with rapid par- 
ticle production - a process known as preheating. During 
preheating, fluctuations of scalar and Gauge fields exhibit 
exponential growth by parametric resonance [^-Q ■ It has 
a host of potentially radical side-effects: Grand Unified 
Scale baryogenesis non-thermal symmetry restora- 
tion |(|, and topological defect formation 0. Here we 
will discuss a side effect which may have persisted until 
the present day - the amplification and sculpting of pri- 
mordial magnetic fields to the amplitudes seen today on 
cosmic scales. 

Magnetic fields are known, partly via the Faraday ro- 
tation of light they induce, to permeate many astro- 
physical systems including intra-cluster gas, quasars, pul- 
sars and spiral galaxies. The fields are large, with mag- 
nitudes ~ 3 x 10~ 6 G on scales greater than 10 kpc M. 



Such amplitudes present an "inverse" fine-tuning prob- 
lem as compared with the standard one in inflation^: 
Since Maxwell's equations are conformally invariant and 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models 
are conformally flat^], the cosmic expansion does not cre- 
ate photons or magnetic fields. The origin of these large 
amplitude fields, correlated on such large scales, is still 
generally regarded as an unsolved mystery, despite the 
proliferation of putative explanations |9|,[l0| . 

The observed magnetic fields today have an energy 
density comparable to that in the Cosmic Microwave 
Background (CMB): r = B 2 /{%tt Pi ) ~ 1 0. If we run 
the cosmic clock backwards past a redshift of z > 100 
where structure formation is strongly in the linear regime, 
r may have decreased to around 10 -34 through the com- 
bined effects of the galactic dynamo |BEo| and collapse of 
structure, which amplifies the magnetic field as (Sp/p) 2 ^ 3 
due to flux conservation. The galactic dynamo efficiently 
converts differential rotation of spiral galaxies into mag- 
netic field energy and without it r ~ 10~ 8 is required to 
seed the observed fields |llj . 

The limit on a homogeneous magnetic field on horizon 
scales today is < 10~ 9 G [T^] . In contrast, at decoupling 
a magnetic field at smaller scales would lead to dissipa- 
tion of energy into the photon fluid and lead to spectral 
distortions. To avoid conflict with COBE FIRAS results 
requires the field to be less than ~ 10~ 8 G today at scales 



^For example, for the potential V = j \<f> 4 , CMB anisotropies 
*With no explicit terms or interactions which break confer- in the absence of preheating demand A ~ 10" 13 , a rather 
mal invariance. severe fine-tuning. 

'i.e., The Weyl tensor, Cap^v, vanishes. 
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0.4 - 600 kpc. 

The time evolution of r is typically believed to be 
rather trivial: r ~ constant. This is due to the high con- 
ductivity of the universe through the matter and radia- 
tion dominated phases which conserves magnetic flux and 



Except for the last section we will use a flat FLRW 
spacetime as a background. The metric is then 



ds 2 = a 2 (i])(-di] 2 + 5 i jdx t dx j ) 



(2.2) 



leads to the behaviour B ~ aT 2 and B 2 / > 7 ~ constant. where V = J a x dt is conformal time, a(rj) is the scale 



However, during preheating and inflation, the low con- 
ductivity of the universe, due to the paucity of charged 
particles, creates an environment in which r can change 
freely. 

The production of magnetic fields during inflation has 
been studied by Turner and Widrow and Davis et 
al. Ilq] and during phase transitions by several authors 
|l4f-jT6| . In reheating their production via stochastic cur- 
rents was investigated by Calzetta et al. 

In this paper we consider the mechanisms discussed by 
Turner and Widrow fllf| and show how preheating may 
lead to resonant amplification of magnetic fields |18|. We 
also discuss a mechanism [|l^,|o| based on the breaking 
of conformal flatness of the background geometry due 
to metric preheating rather than breaking of the confor- 
mal invariance of the Maxwell equations. Although they 
also lead to resonance, we do not consider the axion-like 
couplings 4>F* LV F 11V since they have been considered in 
depth by a number of authors 21 2^]. We will also not 
describe resonant production of magnetic fields in low- 
energy string actions where conformal invariance is bro- 
ken by the existence of the dilaton <j>. Such models have 
been discussed in |2l3 [261]. 



II. 



MAGNETIC FIELDS IN CURVED 
SPACETIME 



Maxwell's equations arise from the Lagrangian density 
— jF^F^u, where F^ v = 2V[ A ,A„] is the Maxwell tensor, 
is the four-potential, is the curved space, covariant 
derivative, and square brackets on indices denote anti- 
symmetrisation on those indices. 

The Maxwell equations that arise are then: 



R llv A v - VuSuA v = , 



(2.1) 



where □ = V^V = {1/^/^)8^ 



-gd u ) and g 



det(g fll/ ). The Ricci tensor term arises through the non- 
commutativity of covariant derivatives and application of 
the contracted Ricci identities 2~S/^A V = R liV A v (27]]. 

The four-potential suffers from a gauge freedom which 
must be eliminated. One may use either the covari- 
ant Lorentz gauge condition VM M = or the com- 
bined Coloumb/tri-dimensional/radiation gauge condi- 
tions A q — 0, d % Ai — 0. In both cases the last term 
in Eq. ( |2.l| ) vanishes |. 



factor of the universe and Sij is the Kronecker delta. The 
traceless part of the Riemann tensor - the Weyl tensor 
Caa^y - defined by WA, 



(2.3) 

vanishes in FLRW backgr ound s which are therefore con- 
formally flat. The metric (2.2) is also conformally static. 

Placing a homogeneous magnetic field in a FLRW back- 
ground is not consistent since the magnetic field picks 
out a preferred direction which is not consistent with 
the maximal symmetry spatial subsections of the FLRW 
models. Instead, the (anisotropic) Bianchi models pro- 
vide an appropriate background for the study of this 
problem [ p8[ . 

Instead we will assume that the magnetic field pro- 
duced will not be coherent on very large scales. Such a 
possibility is already strongly constrained by the CMB. 
Rather we will assume that the power spectrum, B(k), of 
the magnetic field is statistically isotropic and homoge- 
neous, hence consistent with the symmetries of the back- 
ground FLRW model. One then finds, e.g., Q: 

(B z {k)B*{k')) = 4ir 3 5 3 (k - fc')^(fc)|S(fc)| 2 , (2.4) 

where, due to the divB = constraint, P%j{k) must be 
the transverse projection tensor: 



kik ■ 

Pi j W = % p- 



(2.5) 



Assuming the spectrum B(k) is known, then constraints 
at small scales can be used to normalize the spectrum 
and lead to predictions on large scales. 

The energy in the magnetic field in a logarithmic k- 
space interval d In k is 



PB 



d PB _ |B(fc)| 2 



din A; 



(2.6) 



The evolution of magnetic fields is usually described as 
B(k) oc a~ 2 , which means that pb behaves as isotropic 
radiation. 



^If one explicitly breaks the !7(1)em gauge invariance and 



conformal invariance by introducing a photon mass term 
m 2 A V A V into the Lagrangian, then one recovers the Proca 
equation, and the gauge condition V^A^ — 0, becomes a true 
constraint equation. 
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III. A SIMPLE BUT EFFECTIVE ANALYTICAL 
MODEL 

As we shall see, a most efficient and elegant amplifi- 
cation mechanism is to assume a complex scalar field, a, 
charged under U(l), in addition to the inflaton. We will 
assume that its potential, V(aa*), is such that during 
inflation it is displaced from its global minimum. This is 
relatively easy to arrange and occurs rather naturally in 
hybrid models of inflation [|§ 0. 

One way to achieve the desired displacement from the 
global minimum is to give a a negative effective mass 
during inflation which drives it to a non-zero vacuum 
expectation value (vev). At the end of inflation the effec- 
tive mass becomes positive and the field begins coherent 
oscillations. This is a typical scenario for Affleck-Dine 
baryogenesis where the coherent oscillations lead to 
the baryogenesis 

Giving a a non-zero vev during inflation spontaneously 
breaks the U(l) of electromagnetism and causes any 
monopole-anti-monopole pairs to be connected by mag- 
netic flux tubes. These confining flux tubes facilitate 
the annihilation of monopoles. This, the Langacker-Pi 
solution to the monopolc problem, is an independent 
benefit of breaking conformal invariance in this manner. 
Such an additional weapon may be required to deal with 
monopoles produced by non-thermal symmetry restora- 
tion in preheating [|| or in models of inflation which do 
not solve the monopole problem, such as canonical SU (5) 
where the inflaton is a gauge singlet 33 1. 

We will not, however, proceed any further in building 
a detailed phenomenology for a but will assume for ped- 



**For example, consider the archetypal potential |30| 

V — a 2 <j> 2 a 2 + \aaa* — fi 2 \ 2 , 

where a,jj, are constants. Inflation occurs at (f> > </> CI = fi/a 
where the minimum of the potential is \ — X — an d hence 
the effective mass of the photon is zero and the U(l) of elec- 
tromagnetism is unbroken. For — and (f> < (f> CI the 
minimum of the potential now corresponds to the globally 
SUSY vacuum at <f> = 0, x — X* — m/v^- 

"A specific model is given by the following potential in the 
supersymmetric standard model (SSM) along a flat direction 



agogical reasons, to become clear later on, that around 
the global minimum the potential is quartic and the field 
is conformally coupled to the curvature. The Lagrangian 
for this scalar QED is: 



R 



1 



-^{(JCT*) 2 - J^RWl 2 + Anflaton • 



(3.1) 



The conformal coupling will simplify the evolution equa- 
tion for a and reduce it to a form independent of a. The 
gauge covariant derivative = V ^ — ieA^ leads to an ef- 
fective mass for the photon m 2 = 2e 2 |er| 2 which oscillates 
in time as a oscillates. This leads to parametric resonant 
amplification of analogous to studies in Minkowski 
spacetime. 

We work in the so-called unitarity gauge in which a = 
a*, and decompose a into a homogeneous part and a 
fluctuation: a(t, x) — > a(t) + 5<r(t,x). Now let er(ij) be 
the initial amplitude of a oscillations. We assume that 
the oscillations are independent of the inflaton, 0, and 
follow the notation of |3J] in denoting variables rescaled 
by the scale factor a with a tilde; e.g., a = aa. Then the 
equation for a(t) is 



o" + A 



+ e 2 a 2 (A 2 )a = 0, 



(3.2) 



where a prime denotes the derivative with respect to the 
conformal time, 77, and (A 2 ) is the expectation value of 
A^A^. The electromagnetic field vanishes in the back- 
ground, and hence it is automatically gauge invariant in 
the perturbed spacetime. Neglecting the last term in 
Eq. ([T^) for the moment and introducing the dimen- 
sionless quantities x — y/\^a(U)r) and / = a(t)/a(ti), 
we find 



(3.3) 



The solution for this equation can be written as an elliptic 
cosine, / = cn(z, 1/V2), which yields Ull 



(3.4) 



a = a{ti)cn ( x, —= 
v2 



The Fourier modes of a fluctuations satisfy the following 
equation: 



V 



+ h.c. 



+ 



where m a is of order the weak scale, m 3 / 2 is the gravitino mass 
and n is proportional to the number of chiral superfields defin- 
ing the flat direction. During inflation the cH 2 term domi- 
nates and drives a away from the origin. After inflation a os- 
cillates around the time-dependent minimum of the potential. 
The terms proportional to A are soft-supersymmetry-breaking 
corrections responsible for violating B — L and giving rise to 
baryogenesis p2j. 



k 2 + 3cn 2 I x, 



dx 2 

where k 2 = k 2 / (\ a & 2 (t t )) 



1 
71 



5<Tk = . 



(3.5) 



A. Parametric amplification of magnetic fields 



Variations of the Lagrangian ( |3.l| ) with respect to A^ 
leads to the following equation, 
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-f + 2e 2 A v \a\ 2 , 



(3.6) 



where the current f is defined by j v — ie(aVa* — 
ff'W), and vanishes when a = a*. Adopting the 
Coulomb or radiation gauge conditions, Aq = 0, 9M$ = 
0, Fourier modes of A4 satisfy pi|,[l8[ 



(3.7) 



A'l + (k + 2e 2 a 2 a 2 )A k = 



Substituting the solution (J3.4j) for (p.7|), we find 



c?x 2 



2e 2 



A; cn [ X 'T2 



1 



A fc = . 



(3.8) 



0.5 




5 10 15 20 

2e / A, a 

FIG. 1. Density plot of th e Floquet chart for the gen- 
eralized Lame equation (3.8) for < 2e 2 /A CT < 22 and 



K < 3. The shaded regions correspond to parameter ranges 
where parametric amplification of magnetic fields can be ex- 
pected, (if, > 0. The Floquet index, fit, takes larger val- 
ues in the darker shaded regions, and reaches its maxima for 
2e 2 /A ff = 2n 2 at k 2 = 0. 




2e7 A« 



FIG. 2. Three-dimensional plot of the Floquet chart for 
the generalized Lame equation (3.8) for < 2e 2 /A (T < 25 and 
k 2 < 3. 



The whole system reduces to a problem in Minkowski 
spacetime and hence can be sol ved ex actly using the Flo- 
quet theory. In fact Eqs. ( |3.5| ) and (3J3) are the Lame 
and generalized Lame equations respectively. This ele- 
gant exact solution is unstable to perturbations which 
introduce a length scale into the problem (such as giving 
cr a mass) but the existence of the parametric resonance 
is stable. 

The solutions of these equations behave as ~ e^ kX 
where fi k is the Floquet index, which controls the 
strength of the exponential g row th. As for the solutions 
of the dak fluctuation, Eq. (3.5), there is only a single 
resonance band 
Hubble range |35 




constrained to lie in the narrow, sub- 



< k 2 < V3, 



(3.9) 



with a small maximum growth rate, /i max ~ 0.03598 at 
k 2 1.615. In the absence of the a decay to the mag- 
netic field, resonance ends before the energy of the ho- 
mogeneous <7 is sufficiently transferred to the a fluctu- 
ation, in which case the final variance is estimated as 
{5a 2 ) m 0.05a 2 . 



10 
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FIG. 3. The evolution of the variance {A 2 ) = {A^A^} and 
{8a 2 ) for 2e 2 /A CT = 5000 in the Hartree approximation at 
zero conductivity, c.f. Fig. ^| In this case fluctuations of the 
magnetic field are dominated by sub-Hubble modes, and the 
growth of {A 2 ) stops by backreaction effects. Note that {8a 2 ) 
is also amplified with smaller growth rate relative to that of 
{A 2 ). 

In contrast the magnetic fluctuations can exhibit 
strong amplifications, whose strength depends on the ra- 
tio, 2e 2 IXg. According to the analytic investigation in 
Ref. p5[ , the strongest resonance occurs at k 2 = with 
/^max = 0.2377 when the parameter 2e 2 J\ a equals 



2e 2 /A CT = 2n 2 



(3.10) 



where n is an integer. Fluctuations with low momenta 
(k — > 0) are enhanced in the parameter range, 



n(2n - 1) < 2e 2 /A CT < n(2n + 1) . 



(3.11) 
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in which case [ik is typically large and strong resonance 
can be expected. This is found in Fig. [I] where we show 
a density chart of the Floquet index vs k 2 and 2e 2 /A CT . 

When 2e 2 /\ a ~ 0{\) the magnetic field is not 
suppressed on super-Hubble scales during inflation 
S-lil 43 1 . Since the resonance bands (where /ifc > 0) 
stretch down to include arb itrari ly small k/aH in the 
parameter regions given by (3.11), this allows the reso- 
nant production of large-scale, coherent, magnetic fields 
during preheating without violation of causality [ff4|45 
for the case of 1 < 2e 2 /X cr < 3 and 6 < 2e 2 /A (T < 10. 

In Fig. H we plot the evolution of Ak for 2e 2 /X rT = 2 
and a super-Hubble mode k = 10~ 25 . We find that Ak 
is amplified about 10 9 times by parametric resonance, in 
which case the resultant cosmological magnetic field is 
large. However, as we discuss in the next subsection, the 
growth of conductivity during preheating and thermal- 
isation counteracts this resonant growth, and can over- 
whelm it completely. 

For large 2e 2 /A a (~S> 1), the inflationary suppression is 
strong [Q, which makes the large-scale magnetic fields 
negligibly small even if they are amplified by parametric 
resonance. This is actually preferable since development 
of a strong, coherent magnetic field on cosmological scales 
would destroy the isotropy of the background geometry 
set-up during inflation. The magnetic spectrum is blue 
and steep (cx fc 3 ) so that the variance is dominated by 
sub-Hubble modes. 

Since the magnetic field modes are growing expo- 
nentially, backreaction effects become important after 
the fluctuations are sufficiently amplified. Taking this 
into account via the one-loop Hartree approximation, 
Eq. (3.2) is modified to 



+ \ a (cr + 3(Sa 2 })a + e 2 a 2 (A 2 )a = 



(3.12) 



As long as the ratio 2e 2 /A CT lies in the range of Eq. ( [3.11 ) 
the growth rate of A\ is typically larger than that of the 
5(Jk fluctuation.^] This makes to stop the growth of the 
magnetic field modes earlier by ba ckreaction effects when 



the term e 2 a 2 (A 2 )a in Eq. ( 3.1 2| ) is comparable to the 



A<jf7 term, which yields 



(A 2 )^a 2 /(e 2 /\ a ). 



(3.13) 



This relation indicates that the final variance is sup- 
pressed with 2e 2 /A CT being increased, which is similar to 
the standard picture of preheating. In actual numerical 
simulations based on the Hartree approximation, the fi- 
nal vari ance typically takes larger values than estimated 
by Eq. (3.13). The backreaction effect due to the growth 
of magnetic fluctuations do not completely violate the a 



H Note that in the limit of 2e 2 /A CT — > oo, the maximal 
/z*; asymptotically approaches the value 0.2377 for arbitrary 
2e 2 /A ff & 



oscillations p4| , which can lead to the amplification of 
the a fluctuations even after magnetic fluctuations are 
sufficiently amplified. This behaviour is found in Fig. || 
where we plot the evolution of fluctuations for the case 
of 2e 2 /A ff = 5000. 



B. The growth of conductivity, <j c 

The above analysis assumed that the conductivity a c 
of the universe vanished during inflation and preheating. 
This is almost certainly incorrect |3f| but accurate mod- 
elling of the growth of conductivity is difficult for two 
reasons: 

(1) such calculations depend sensitively on the under- 
lying theory in which the inflaton is embedded, and 

(2) to estimate the growth of conductivity a c requires 
non-perturbative, non-equilibrium quantum field theory 
techniques, hence is extremely difficult. 

(3) Accurate estimates of the final magnetic field re- 
quires the conductivity in three phases - during inflation, 
during the initial resonance phase and during thermali- 
sation, each of which is dominated by different physics. 

The growth of a c in the QED case has been studied 
in detail p9j . Since this is not appropriate for energies 
near the GUT scale we can only draw broad lessons: the 
conductivity grows exponentially but is also spatially in- 
homogeneous. 




(I 50 100 ISO 200 

X 

FIG. 4. The growth of Ak in the case of 2e 2 /A CT = 2 for 
the mode k — 10~ 25 . In the absence of conductivity (a c = 0) 
the large-scale magnetic field fluctuations are strongly am- 
plified due to the oscillating a field. As the rapidity of the 
growth of the conductivity increases (increasing c) the reso- 
nant growth of At is increasingly stalled. For large conduc- 
tivity (not shown) a c — > oo the electric field vanishes and the 
magnetic field becomes frozen, Bk oc a~ 2 . 

This is related to the fact that while the plasma is on 
average charge neutral, there will be fluctuations in the 
charge density which act as stochastic sources of mag- 
netic fields; see p7| . Given the problems described above 
we take a phenomenological approach to the growth of 
conductivity. 
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Since we are interested in large scales we neglect the 
spatial variation of a c and model its growth as: 



o c = <7f tanh[c(a; — Xq)\, 



(3.14) 



where <jf is the final value of conductivity, and c controls 
the growth rate of a c , x is the dimensionless conformal 
time and xq is the onset of preheating. We therefore 
assume that a c — during inflation. As noted in |38| if 
<7 C > e 2 er 2 at the onset of preheating then the resonance 
in Ak never begins. 

a f and c determine the strength of conductivity. Fig. |] 
shows the evolution of a cosmological Ak mode for three 
pairs (tJ/,c). The value cf/ = af /(^/X^a(ti)) = 1 is used, 
where a(ti) is the value of the a field at the beginning of 
its oscillations. 

The four-potential in the finite-cr c case obeys the equa- 
tion: 



(k 2 + 2e 2 a 2 a 2 )A k = -a c aA' k , 



(3.15) 



which shows how the conductivity acts to damp the res- 
onance when we neglect the spatial dependence of a c . 
Fig. ^ shows how the the preheating resonance competes 
with the damping due to conductivity. If the growth of a c 
is too rapid (large c) the resonant growth of A k is stalled. 
However, for relatively slow growth of <j c (roughly more 
than 10 a oscillations) A k can grow almost to its a c = 
maximum. 

Once backreaction causes the resonance to end there is 
nothing to compensate the damping effects of the finite 
conductivity and A k begins to decay expo nentially. In 
the limit a c — > oo the solution of Eq. (3.15) is obviously 
A k = constant, which corresponds to B k oc a~ 2 , the ratio 
of magnetic field energy density to incoherent radiation 
energy density is fixed. 

To estimate r at the start of the radiation dominated 
phase is therefore a subtle issue because one must know, 
not only the growth of conductivity during the initial 
preheating phase, but also how the conductivity grows 
during thermalisation and how a decays by Born process 
to complete reheating. If the conductivity is high during 
preheating, the magnetic fields will exhibit exponential 
suppression during which a c increases from zero to the 
final value, cry, which means that the gains of preheating 
will be washed out and lost. 

When the conductivi ty ter m is much smaller than the 
2e 2 a 2 a 2 Ak term in Eq. ( [3.11: ) during preheating, the evo- 
lution of magnetic fields is the same as the case of the 
non-conduc tivity in preheating phase. However when the 
rhs of Eq. ( |3 .15 ) becomes of order the 2e 2 a 2 a 2 A k term 
after preheating, the magnetic field begins to be expo- 
nentially suppressed. 

Although the A k freeze when the lhs of Eq. ( [3.15 ) be- 
comes negligible relative to the conductivity term, the 
gains obtained in preheating are not generally preserved 
due to the rapid decay of magnetic fields before the freeze 
of Ak- However, the Born decay of a before thermalisa- 
tion can alter the strength of the 2e 2 a 2 a 2 Ak term, which 



may alter the above estimates. In addition to this, we 
need to know the evolution of conductivity during ther- 
malization for a complete study, although it is difficult 
and few studies of thermalization after preheating exist; 
see e.g. M. 




FIG. 5. The maximum variance {A 2 } as a funct ion o f the 
conductivity parameters c and 5f appearing in Eq. (3.14). As 
c and bf increase the period and degree to which preheating 
can amplify (A 2 ) decreases rapidly. 

In conclusion, in the absence of conductivity, either the 
magnetic field is resonantly amplified on super-Hubble 
scales, or it has a k 3 spectrum and is too small on cos- 
mological scales. 

When conductivity is included one introduces several 
model-dependent parameters into the problem which im- 
pact on the viability of preheating as a significant source 
of magnetic fields. The effect of resonance is typically 
washed out by the growth of conductivity, while the 
final size of magnetic fields depends on details of the 
Born decay process and the evolution of conductivity af- 
ter preheating. To answer which case applies is model- 
dependent and beyond the scope of the current paper. 



IV. GEOMETRIC MAGNETISATION 

Turner and Widrow jll[] found that the most efficient 
way to produce magnetic fields is to break U(l) gauge 
invariance as well as conformal invariance, via the La- 
grangian: 



L = 



R 



16ttG 

7 



\f^ v - \ra^ 



2 R^A»A* 



inflaton ; 



(4.1) 



where /3, 7 are real constants. These terms were found 
to give rise, for apparently reasonable values of the con- 
stants, during inflation, to fields corresponding to the 
required value r ~ 10~ 8 today Variation of the ac- 
tion (4T) with respect to A^ and P" yield the following 
equations of motion: 



V% - 0RA u - 7 fl£4u = 0, 



(4.2) 



G 



9A + d x F, y + d v F; 



\f.i 



0. 



(4.3) 



Writing these equations in terms of magnetic and electric 
fields and eliminating the electric field, we obtain 



with 



(a 2 B)" - V 2 (a 2 B) + %)a 2 B = 0, 



a a V a 



(4.4) 



(4.5) 



Expanding the magnetic field in Fourier components as 
a 2 B = J e _tk x Kfc G? 3 k, each mode satisfies the following 
equation, 



+ \k 2 + 



9{rj)] B k = 0. 



(4.6) 



Hereafter we set 7 = and leave (3 as a free parameter, in 
which case 9 reduces to 6 = 6(3a" ja = &(3a 2 R. When the 
system is dominated by the inflaton field, 0, the scalar 
curvature is: 



R = 



8tt 



'pi 



(4.7) 



where V((f>) is the inflaton potential. During inflation R 
slowly decreases. When (3 is negative, the magnetic field 
fluctuations exhibit super-adiabatic amplification due to 
the so-called negative coupling instability, as studied in 
the non- minimally coupled case in Refs. This 
enhancement is most relevant during inflation. 



B 



k oca 



(4.11) 



for negative values of (3. In this case the energy density 
in the /c-th mode of the magnetic field evolves as 



PBO c\B k \ 2 /a 4 = aV r - 



48/3-5 



(4.12) 



which means that the ratio, r = ps/pj, increases during 
inflation when (3 < 0. This makes it possible to reach 
the value r ~ 10 -8 required to explain the existence of 
current galactic magnetic fields pJ§ . 

Large negative values of (3 lead to extremely strong am- 
plification of magnetic fields. When [3 = —1/6, B k oc a 
and pb oc a~ 2 , which corresponds to the minimally cou- 
pled scalar field case. Compared with the standard adi- 
abatic result, ps oc a~ 4 with (3 — 0, the energy density 
decreases more slowly due to superadiabatic amplifica- 
tion. 

For (3 < 1, super-Hubble magnetic fluctuations ex- 
hibit enormous amplification during inflation, i,e., r — 
Pb/ oc a c with c > 6, which conflicts with observations 
unless their initial values at the start of inflation were 
extraordinarily small. 

When (3 is positive, magnetic fields are exponentially 
suppressed during inflation. For < (3 < 1/48, which 
corresponds to < v < 1/2, B k and pb evolve as 



Sfcoca-^V 1 ^- 1 ), pBoza-V^^s. (4.13) 
When (3 > 1/48 (i.e., complex u), we find 



B k ex a 1 ^ 2 , oc a 



(4.14) 



A. Magnetic amplification during inflation 

Let us study the evolution of magnetic fluctuat ion s 
during inflation. When 7 = 0, the solutions for Eq. ( |4.6| ) 
are expressed as combinations of the Hankel functions 
1,2) H : 



Hl J) (J 



B k = dy/rjH^ikri) + c 2 ^H^\kri), 



(4.8) 



where c\ and C2 are constants, and the order v of the 
Hankel functions is given by 



v 2 = - - 12(3. 



(4.9) 



The choice of c\ — and c-i = corresponds to 

the Bunch-Davies vacuum. In the long wavelength limit, 

krj — * 0, HI, ' (krj) approaches the values 



(4.10) 



where F(^) is the Gamma function. In inflation, confor- 
mal time can approximately be written as rj pz —l/(aH), 
and long wave B k modes exhibit exponential growth, 



in which case the evolution of magnetic fields is indepen- 
dent of the strength of (3. 
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mt 

FIG. 6. The evolution of a super-Hubble Bk mode during 
inflation and preheating for (3 = 0, —0.01, 10 in the massive 
chaotic inflationary model. We choose the initial value of in- 
flaton as 4> = 3m p i, which corresponds to 60 e-foldings before 
the end of inflation. The inflationary period continues until 
mt ~ 20, after which the system enters the reheating stage. 
While Bk is exponentially suppressed during inflation for pos- 
itive (3, negative (3 leads to superadiabatic amplification. 



7 



We plot in Fig. ^ the evolution of a super-Hubble Bk 
mode for (3 = 0, —0.01, 10 for the inflaton potential, 



v{4>) = z m 



2,2 



(4.15) 



Whe n (3 — 0, Bk is constant (i.e., Bk oc a ) from 
Eq. ( 4.13| ), as is confirmed in Fig. |[ For negative 
(3, B k exh ibits an exponential increase as estimated by 
Eq. (4.11). An important point to note is that the rapid 
growth of magnetic fields may affect the evolution of 
background quantities, an effect we do not include. 

In Ref. p8 |, it was found that exponential growth of 
scalar field fluctuations makes the inflationary period ter- 
minate earlier, in the context of a non-minimally cou- 
pled scalar field. Exponential growth of large scale mag- 
netic fields will also stimulate the enhancement of super- 
Hubble metric perturbations, which may lead to devia- 
tions from the scale-invariant Harrison-ZePdovich spec- 
tra. This is expected to be strong for (3 < — 1 in the 
analogy of the non-minimally coupled scalar field case 
p8| . Complete analysis including backreaction and met- 
ric perturbations is now in progress. 



B. The preheating phase 

The inflationary period corresponds to < mt < 20 in 
Fig. [| after which time the system enters the reheating 
stage. During reheating, the scale factor ev olves asooc 



cx rf in the massive inflaton potential (4.15). From 



t 2 / 3 , . 

Eqs. and (pLlcb , we have that B k cx a3 (1+ V 1 " 48 ' 3 ) 

and pb oc a'i^ 1-48 ' 3-7 -' for negative (3 in the long- wave 
limit ki] — > 0. Similarly Bk oc a 1 / 4 and ps oc a -7 / 2 when 
(3 > 1/48. However, this corresponds to an estimate 
of the frequency 8 oc r\~ 2 oc £~ 2 / 3 , which only provides 
information about the average amplitude of the scalar 
curvature. 

In actual fact the scalar curvature oscillates due to the 
oscillating inflaton condensate, which can lead to effi- 
cient enhancement of field fluctuations 49 - 5j] . We find 
in Fig. H that Bk begins to grow for mt > 20 in the 
case of (3 — 10 in spite of the inflationary suppression. 
This is the geometric preheating stage where Bk grows 
quasi-exponentially, in which case the above naive esti- 
mate neglecting the oscillations of the scalar curvature 
can not be applied. In the non-minimally coupled multi- 
field case, the growth of scalar field fluctuations during 
preheating is only relevant for |£| > 1 p9|,p0|. 



FIG. 7. See associated JPEG image. The 

time-averaged Floquet chart for the geometric magnetisa- 
tion mechanism during preheating. Since the problem is not 
scalar- factor independent, exact Floquet theory cannot be 
used - the expansion causes the weakening of the resonance 
bands but also removes the stability bands so that all modes 
are amplified for sufficiently large (3. It is clear that as (3 is 
increased, the super-Hubble modes k < 1 are the ones that 
grow first. The negative coupling case (3 < 0, which is much 
stronger, is not shown. 

Let us analytically study the evolution of magnetic 
field fluctuations during preheating. Making use of the 
time averaged relation, (t;<() 2 )t = (V(4>))t, with the po- 



tential (4.15), the evolution of the inflaton condensate is 
described by 



m p i 



Sirmt 



■ sin mt, 



(4.16) 



where we choose the time when the oscillation starts as 
mf o = 1/4 || . Then the scalar curvature fl4.7| ) is approx- 
imately 



R 



M 2 



(1 — 3 cos2mi). 



(4.17) 



Although R oscillates, its amplitude decreases as t~ 2 due 
to the cosmic expansion which means that parametric 
resonance soon becomes ineffective if \(3\ is small. Sub- 



stituting Eq. (4.17) into Eq. (4.6) and introducing a new 
scalar field Bk — a x l 2 Bk, Bk satisfies the well known 
Mathieu equationPl 



d 2 B k 
dz 2 



{A k -2qcos2z)B h =0, 



where for (3 > 0, 



A k = 
and for (3 < 0, 

Ai, = - 



3 q 



3(3 

tt 2 * 2 ' 



;1 



q 



2+2 ' 



7T t 



(4.18) 



(4.19) 



(4.20) 



Here z = mt, and t — mt/ '(2tt) naively corresponds to the 
number of oscillations executed by the inflaton at time t. 

In the context of standard preheating with the effective 
potential, V(<j),x) = \m 2 § 2 + ^g 2 <ft 2 x 2 , the relation of 
Ak and q for the x ne ld is written as Ak = 2q+k 2 / (m 2 a 2 ) 
Q . In this case x particle production is inefficient unless 
the initial q is much larger than unity. In contrast, the 
resonance band is broader in the present model BOf . 



^We neglect the ( a/2a ) 2 — a/2a term which appears in the 
parenthesis of Eq. (4.18), which can be justified for |/3| > 1. 
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Therefore parametric resonance takes place for smaller 
initial values of q. In spite of this, since the resonance 
band is narrow for q < 1, we typically require the cou- 

(see Fig. 



> 1, for relevant growth of Bk 



pling, 

where we show Floquet indices for positive (3). When 
|/3 1 > 1, fluctuations grow as Bk ~ e Mfemt , whose growth 
rate gets gradually larger with \f3\. For \(3\ 3> 1, however, 
the final variance of magnetic fields will be suppressed as 
studied in e.g., Ref. ]50| ]. 

When (3 is positive, long-wave Bk modes are exponen- 
tially suppressed during inflation. Hence it is rather diffi- 
cult to produce sufficient large scale magnetic fields even 
when (3 ^> 1. On sub-Hubble scales, however, the in- 
flationary suppression of Bk is weak relative to super- 
Hubble modes and magnetic fields are excited during 
preheating. Hence the final magnetic variance (A 2 ), is 
dominated by sub-Hubble modes. For negative values of 
(3 with (3 < — 1, the growth of Bk can be strong but 
is typically dominated by the growth during inflation. 
When — 1 < (3 < 0, the production of magnetic fields is 
weak during preheating, while they are amplified in the 
preceding inflationary phase as found in Fig. |6|. 



C. Effects of the growth of conductivity 

Now let us consider the ratio r = Pb/p^ on some co- 
moving scale in the presence of the preheating phase. As 
the reheating process proceeds, the effect of the conduct- 
ing plasma is expected to be important jy]. This effect 
appears as a friction-like term in the equation of Bk 



BZ+[k 2 + 8( V )] Bk = -a c aB'k 



(4.21) 



where a c is the conductivity of the plasma. If the conduc- 
tivity is very high, we find Bk ~ constant, which implies 
that the energy density of magnetic fields decreases as 

-4 

Pb ~ a 

We assume that the effect of the conductivity begins to 
dominate at some temperature, T c (< T r ), where T r is the 
reheating temperature. At the first Hubble-crossing dur- 
ing inflation, the ratio of ps to the total energy density, 
Pt, is approximately estimated as pb/pt ~ {M/m p \) A , 
where M 4 is the energy scale of inflation. For negative 
(3, one obtains the following ratio on the comoving length 
scale A neglecting the parametric amplification of mag- 
netic fields during preheating 1 1 1 : 



10 *-., ( ".)""( JL) 

\m pl J \m p ij 



(l+3ot)/3 



2k 



-4(l+a)/3 



( X 



VMpc 



a — 5 



(4.22) 



where a = yT — 48/3. For example, when (3 = —1/2, 
r|iM P c reaches - 1CT 8 for M = 10 17 GeV, T r = 10 17 
GeV, and T c = 10 17 GeV, in which case seed magnetic 
fields can be produced without the need for the galactic 
dynamo mechanism. When — 1 < (3 < 0, since parametric 
excitation of magnetic fields is irrelevant, the estimation 
of r in Eq. ( 4.22Q is hardly modified due to the existence 
of the preheating phase. 

In contrast, for (3 < — 1, it is expected that preheating 
will lead to the increase of r. In this case, however, r|iM pc 
is typically much greater than unity even in the abse nce 
of preheating because lO 26 ^" 5 ) > 1 in Eq. (fO|f^ 



Although r is further increased corresponding to the am- 
plification of Bk during preheating, this case will be ruled 
out by observations. 

Let us consider the case j3 > 1 where the excitation of 
magnetic fields by resonance is expected. In this case, an 
analytic estimate of r neglecting the contribution during 
preheating is 



10 



-130 



(M_\ ( Tr\ 



1/3 



\m p \J \m p \J 



(—) 



-4/3 



A 



Mpc 



(4.23) 



Due to the strong inflationary suppression, r is restricted 
to be very small. For example, for M = 10 16 GeV, 
T r = 10 15 GeV, and T c = 10 15 GeV, r|i Mpc = 1(T 129 . 
During preheating, the Bk fluctuation exhibits exponen- 
tial i ncrea se, which makes r larger than estimated in 
Eq. ( 4.23 ). For example, when (3 — 10, Bk is ampli- 
fied about 10 5 times (see Fig. ||), and the ratio increases 
to r|iMpc ~ 10~ 120 . However, the amplification during 
preheating in the positive f3 case is typically insufficient 
to explain the large-scale seed magnetic fields even for 
/3> 1. 

We conclude that with regard to the geometric mag- 
netisation mechanism, the ratio r = pb/pt is mainly 
determined by the inflationary phase, despite the fact 
that magnetic fields can be amplified during preheating. 
While we have studied this in the massive inflaton model, 
we expect similar results in other inflationary models. 
For example, in the quartic inflaton potential the fre- 
quency 9 depends explicitly on the scale factor and we 
cannot reduce the problem to one in Minkowski space, as 
we did in section III. 



m For example, for /3 = — 1, M 
and T c = 10 15 GeV, we find r 
excessive. 



10 16 GeV, T r = 10 15 GeV, 



IMpc 



10 37 , which is clearly 



***When \(3\ = 1, the initial value of q is estimated to be 
qi w 4.9 with ii = 1/4. 







V. MAGNETIC FIELD AMPLIFICATION DUE 
TO LARGE METRIC PERTURBATIONS 

Since the FLRW metric is conformally flat, i.e., the 
Weyl tensor vanishes, magnetic fields are not produced 
due to the cosmic expansion. During preheating however, 
scalar metric perturbations can grow exponentially on 
both super-Hubble and sub-Hubble scales j34],[l5| . 

This growth of metric perturbations means that the 
spacetime may no-longer be well-described by a confor- 
mally flat background metric. If the metric perturba- 
tions remain small, this breaking of conformal invari- 
ance is small (as measured by the curvature invariant 
Ca/3/ii/C Q/3A " y ) and the production of photons is very sup- 
pressed. Once the metric perturbations at a certain scale 
become large, however, this is no-longer true and the 
production of magnetic fields can be expected. This was 
discussed by Calzetta and Kandus (5i| in the context of 
structure formation and suggested in the context of pre- 
heating in Here we follow the recent analysis of 
Maroto §0|. 

The line element for a flat FLRW model with scalar 
metric perturbations in the conformal Newtonian or lon- 
gitudinal gauge is |5^|34j 

ds 2 = a 2 (r))[-(l + 2<J>)d?7 2 + (1 - 2$)6 ij dx i dx j ]. (5.1) 

We consider the following two-field model in the presence 
of magnetic fields: 



L = 



16ttG 4 ^ 2 V v > 4 v 



(5.2) 



where x is a scalar field coupled to inflaton, (j). Then the 
magnetic field satisfies the Maxwell equation, V — 
0, i.e., 



8 
82^ 



-99» a 9 u0 id a Ap - dpA a )] = 0. 



(5.3) 



Using the relations yf—g = a 4 (l — 2$), g 00 
2$), of = a~ 2 (l + 2$) in the perturbed metric ( |5.l| ) 



-a- 2 (l 



Eq. (5.3) yields for v = i 
8 



dr) 
8 



[(1-2$) (diAo-doAi)} 



+ _[(l + 2$)(8 i A i - 8^)1 = 0. 



(5.4) 



Adopting the Coulomb gauge condition: A = 0, d % A i = 
0, one finds that 



A 



W 2 Ai = 2<f>'A'i + 4$V 2 A l 
8$ (dA t dA 3 



(5.5) 



The effect of metric p ertu rbations appears at second or- 
der in the rhs of Eq. (5.5). In Fourier space this leads to 



convolutions of the form J d 3 k' $ k , A k _ k , , which lead to 
mode-mode coupling. 

However, if we assume that $ is only dependent on 
time on scales larger than some cosmological scale A c = 
2ir/k c ^0|, each Fourier component of Ai satisfies the 
simple equation: 



k 2 A k = 2&A' h - 4fc 2 $A 



(5.6) 



where the coupling between the metric potential on 
smaller scales (k > k c ) and the magnetic fields are ig- 
nored. Note that one cannot simply assume $ = <E>(i) on 
all scales since then the Weyl tensor vanishes identically 
and no photons are produced. 

Treating the full problem is compl icat ed due to the fact 
that the last term in the rhs of Eq. (5.5) does not vanish, 
hence the various components of Aj are coupled. While 
the precise analysis including these fully nonlinear effects 
is very complicated, we can still estimate the amplitude 
of magnetic fields produced during preheating by using 
Eq. (Q. 

Introducing a new fiel d, Ak = (1 — &)A k , to eliminate 
the A' k term in Eq. (5.6), one finds 



k 2 A k = 



(5.7) 



where we have neglected the last term in Eq. (|5.6| ) which 
is not important on large scales. Before the start of pre- 
heating, the A k term is negligible and A k is described 
by the following positive-frequency solution: 



1 



,-ikri 



2k 



(5.8) 



One finds the solution for Eq. (5.7) in integral form |55|: 
A h {r,) = A® + - [ V $"A fe (r/) sin k( v - rf)drf . (5.9) 



The energy density in the magnetic field can be expressed 
as 

PB = {k/af\p k \ 2 , (5.10) 
where the Bogolyubov coefficients, f3 k , are approximately 



Pk = -i 



(5.11) 



Substituting Eq. ($>S$ for Eq. (|5.1l[) with Eq. (jSJf and 



assuming that $' vanishes before and after preheating 
[i.e., &(ti) = and = where the subscript i 

and / denote the values at the beginning and end of 
preheating, respectively], one easily finds that the next 
order term in Eq. (5.9) gives an important contribution 
to j3 k , yielding M 



k 



i 

2k 



fit 



(5.12) 
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where we considered th e sup er-Hu bble m odes: « 1. 

fla.lCj ), and ( 5.12| ), we obtain the 



Combining Eqs. (2.6) 
amplitude of magnetic fields as 



\B k \ « * / (Vfdv- 

Jrn 



(5.13) 



In order to analyze the evolution of the magnetic fields, 
it is convenient to rewrite Eq. ( |5.13 ) using the dimen- 
sionless conformal time, x = y/Xaifar), as 



\Bi. 



(5.14) 



At the decoupling epoch where the coherence scale corre- 
sponds to (k/ ^)dcc ^ 10 ^ GeV, the amplitude of mag- 
netic fields can be estimated by 



Odec J Xi \dx J 



(5.15) 



where we used the value, VX(f>i ~ 10 12 GeV. The ra- 
tio ai/adec depends on the reheating temperature, Tr. 
If the energy of inflaton at the end of inflation were 
instantaneously transferred to radiation, the reheating 
temperature would be Tr ~ 10 15 GeV, which yields 
ai/a<iec ~ Tdec/Tii ~ 10~ 25 . Note that the ratio ai/a^cc 
becomes larger for lower reheating temperature. 

Primordial seed magnetic fields for the galactic dy- 
namo mechanism are in the regions of |-B^ CC | = 10 -25 
G ~ 10~ 15 G. In the single field case in which large- 
scale metric perturbations are hardly amplified during 
preheat ing, i t was found that magnetic fields estimated 
by Eq. ( 5.15 ) are below the values required for the galac- 
tic dynamo in the realistic values of a, / adoc p0| ] ■ 

In the two-field case with self-coupling inflaton, we 
can expect the growth of metric perturbations due to 
the enhancement of field perturbations, which stimu- 
lates the growth of magnetic fluctuations through grav- 
itational scattering. Decomposing the scalar fields as 
ifj(t,yi) — > <fj(t) + S(pj(t,ic), the Fourier transformed, 
perturbed Einstein equations are 



3H6<p k + 
4# fc + 21 



- + 3X(<p 2 + (5^))+g 2 ( X 2 + (6 X 2 )) 



growth of <&k and 8<j) k on la rge s cales, as is found in nu- 
merical simulations of Eqs. ( 5.18 ) and ( 5.16 ). Neglecting 
metric perturbations on the rhs of Eqs. ( 5.17| ) which are 
small during inflation, we find the following analytic so- 
lution 



with 1 34 1 



2^2 



" =4" 



H 2 



2ttX 



(5.19) 



(5.20) 



since H 2 « 27rA</> 4 /3. In the centre of the first resonance 



band, g /X = 2, v is negative only when 4> < 2/vo7r ~ 
0.7m p i, which means that the exponential suppression 
can be avoided during most of inflation. In this case 
large-scale metric perturbations are significantly ampli- 
fied during the preheating phase. 

In Fig. |8| we plot the evolution of $> k , S Xk , ^Xfc: an d 
M(x) = J {d<&/dx) 2 dx for g 2 /X = 2 during inflation 
and preheating for a cosmological mode. We include sec- 
ond order field and metric backreaction effects as spatial 
averages for background equations (see Refs. |43| , [5|j56[ 
for details), and choose initial values for the scalar fields 
at the start of inflation to be 0(0) = 4m p i and x(0) = 
10~ 3 m p i. Metric perturbations begin to grow during pre- 
heating after S Xk grows to or order 5<fi k , which results in 
the final amplitude of order $ ~ 0.1, clearly in conflict 
with observations of the CMB. 

In spite of this, it is worth investigating this case in 
order to understand how the growth of metric pertur- 
bations affects the evolution of magnetic fields. The 
M(x f ) = II s (d$/dxf dx term on the rhs of Eq. fl5.15| ) 
becomes of order 0.01 (see Fig. ||), and the resulting 
magnetic field at decoupling is then estimated to be 
\ B k e °\/ 1G ~ ICTWadcc- When a,i/a dec > 10~ 21 which 
corresponds to the reheating temperature, Tr < 10 11 

GeV, magnetic fields exceed the value, \B k icc \ — 10~ 25 
G, which is required to seed the galactic dynamo. 



S Xk + 3B5xk + 



3H^ k ~2g 2 ^ X S Xk , 



+ g 2 ( C f + (5<j ) 2 )) 8 Xk 



8<pk 

(5.16) 



A X <$> k + 2{ X + ?,H X )<S> k -2g< 



$ fc + H$ k = 4irG(<p6(j> k + X Sxk) 



(5.17) 



(5.18) 



As long as S Xk fluctuations in low momentum modes are 
not strongly suppressed during inflation (i.e., g 2 /X < 10) 
and are excited during preheating, this leads to the 
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FIG. 8. The time evolution of 8\k, 8(j>k, and 

M(x) = (dQ/dx) 2 dx for g 2 /X = 2 during inflation and 
preheating for a cosmological mode. In this case, the en- 
hancement of metric perturbations leads to the production of 
magnetic fields due to the breaking of conformal flatness of 
the background metric. 



With the increase of g 2 /X, the inflationary suppres- 
sion for long wave Sxk modes begin to be signifi- 
cant. For example, in the centre of the second reso- 
nance band, g 2 j\ = 8, the suppression is relevant for 
<f> < 1.3m p i. In the Hartree approximation, the en- 
hancement of super-Hubble metric perturbations during 
preheating was found to be weak for g 2 /X > 8 due to 
the suppressed \ fluctuation at the end of inflation Q , 
which means that magnetic field fluctuations are hardly 
enhanced by large-scale metric perturbations. 

However, since sub-Hubble 5\k fluctuations are free 
from strong inflationary suppression and exhibit para- 
metric amplification during preheating, metric preheat- 
ing is typically vital on sub-Hubble scales j56|. Then 
the mode-mode coupling between sma ll-scale metric and 
large-scale magnetic field in Eq. (5.5) may lead to the 
production of magnetic fluctu ations. In this case ana- 
lytic estimations by Eq. ( |5.13 ) can no longer be applied, 
and we have to solve the complicated nonlinear equation 
(5.5) directly. Whether magnetic fields can be sufficiently 
amplified by the growth of small-scale metric perturba- 
tions is uncertain at present. We leave to future work for 
the precise analysis of this issue. 

We should also note that parametric excitation of sub- 
Hubble Sxk modes will stimulate the growth of large-scale 
54>k and modes. The Hartree approximation misses 
this rescattering effect |37],[38|, which is expected to be 
important once fluctuations begin to be amplified sig- 
nificantly. In fact it was recently found that rescattering 
can lead to the amplification of super-Hubble metric per- 
turbations even for g 2 /A > 8 in one-dimensional lattice 
simulations fl59| . It is unknown whether this holds true 
for g 2 /\ ^> 1, which will be clarified by fully nonlinear 
three-dimensional calculations. 

It is certainly of interest to find parameter regions 
which satisfy both the CMB constraints and produce 
sufficient large-scale seed magnetic fields. Although we 
have restricted ourselves in the chaotic inflationary sce- 
nario, the ratio ai/adec and the energy scale of infla- 
tion are model-dependent. It is encouraging that we 
can test inflationary models by the magnetic fields pro- 
duced, together with CMB and primordial black hole 
over-production constraints during preheating frM]5(i] . 



VI. CONCLUSIONS 

In this paper we have considered the amplification 
of (hyper-)magnetic fields during inflation and preheat- 
ing. The conformal invariance of the standard Maxwell 



equations and the conformal flatness of the FLRW back- 
ground leave the observed cosmic magnetic fields as a 
major mystery. In order to overcome such obstacles, we 
have considered three specific mechanisms: 

(1) Couple the magnetic field to a coherently oscil- 
lating scalar field which induces resonant growth of the 
magnetic field. In the presence of plasma effects, para- 
metric amplification of magnetic fields is typically coun- 
teracted by the growth of conductivity. This competition 
is model dependent and the final outcome depends sensi- 
tively on the conductivity during inflation, the resonance 
and thermalisation phases (see Figs 4,5). 

(2) Break conformal invariance of Maxwell's equations 
through non-renormalisable couplings to the curvature 
such as RA^A^ . When the corresponding coupling con- 
stant, (3, is negative, strong amplification of the magnetic 
field occurs during inflation. As a result it is a promis- 
ing mechanism, though some fine-tuning may be required 
not to over-produce the magnetic fields by the end of pre- 
heating. For positive (3 the produced field is too weak to 
be relevant even with the resonance from preheating. 

(3) Break the conformal flatness of the background 
metric. During metric preheating super-Hubble metric 
perturbations grow exponentially. The resulting growth 
of the Weyl tensor leads to amplification of the magnetic 
field, which while it is generic, is a complex, mode-mode, 
coupling problem. 

It is certainly of interest to consider issues such as the 
non-equilibrium aspects of the problem and a detailed 
model of e.g., the GUT gauge group and couplings be- 
tween the relevant gauge fields and the curvature/other 
fields, which we leave to future work. 
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APPENDIX: MAGNETIC FIELDS WITH RF^F^ 
INTERACTIONS 

The 1-loop QED result in curved space includes terms 
of the form RF^F^ together with similar terms involv- 
ing R^ and R^ a f3- These are more complex to treat as 
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resonance systems because of periodic divergences. To 
illustrate this we consider the Lagrangian 

C = — ^ - ~ ( 1 + b-^) F^F^ + A nfla to„, (6.1) 

167rG 4 \ m^J 
where & is a constant P 



The equation of motion for the Fourier modes of 
are 



bR lk 



0. 



(6.2) 



In the limit of R m 2 (the one appropriate for the early 
universe Jill 
is independent of b 



the coefficient of A' ik becomes R' /R and 



Since R oscillates through zero [see e.g., Eq. (4.17)], the 
equation is not amenable to simple numerical analysis. In 
this regard it is similar to the evolution equation for the 
potential <& in the single, oscillating, scalar field case 0. 
As discussed at the end of Q], the periodic singularities 
do not forbid resonance bands. In the case of negative 
b the possibility of efficient amplification during inflation 
exists due to the negative coupling instability. 
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